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Abstract 



We study the magnetic properties of the coloured multiquark states in the quark- 
gluon plasma where the gluons are deconfined and the chiral symmetry is still broken, 
using the Sakai-Sugimoto model. There are two possible magnetized multiquark con- 
figurations. Both configurations converge to the same configuration at the critical 
field and temperature before they dissociate altogether either into less coloured multi- 
quarks or into other phases for a fixed density. It is also found that the multiquarks 
with higher colour charges respond more to the external magnetic field in both the 
magnetization and the degree of chiral symmetry breaking. Magnetic field also makes 
it more difficult for multiquark states with large colour charges to satisfy the equilib- 
rium condition of the configuration in the gravity dual picture. As long as the chemical 
potential > fionset , the magnetized multiquarks phase is thermodynamically preferred 
over the magnetized vacuum. Pure pion gradient and the chiral-symmetric quark-gluon 
plasma (xs-QGP) phase for the general Sakai-Sugimoto model are discussed and com- 
pared with the multiquark phase in the presence of the magnetic field. It is found that 
at large densities and moderate fields, the mixed phase of multiquarks and the pion 
gradient is thermodynamically preferred over the xs-QGP. 
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1 Introduction 



There has been increasing interest in the study of nuclear phase structure as well as properties 
of a number of nuclear phases, especially the quark-gluon plasma in the recent few years. 
This is due to the new perspective in the nature of strongly interacting gauge theory from 
the holographic principle. Motivated by the AdS/CFT correspondence [H E], a number 
of gravity models was constructed to provide shadow gauge theories which share certain 
essential features with the QCD in the strong coupling regime. Sakai and Sugimoto [3l H] 
proposed a toy holographic model of QCD where chiral symmetry breaking can be addressed. 
In Sakai-Sugimoto model, gluon deconfinement and chiral symmetry restoration are two 
distinct phase transitions. For non-antipodal case, the chiral symmetry restoration occurs 
at higher temperature than the gluon deconfinement [S], therefore it is possible to have a 
nuclear phase where gluons are deconfined while the quarks and antiquarks could still form 
colour bound states. 

Bergman, Lifschytz, and Lippert [6] shows that when the baryon density is sufficiently 
large and the temperature is not too high, gluon-deconfined phase with broken chiral symme- 
try accommodates a nuclear phase where baryons can exist with thermodynamical stability. 
Even though the baryons can exist within the phase, the quark matters containing only 
free quarks or antiquarks do not share the same thermodynamical stability. This can be 
understood as a sign of chiral symmetry breaking, the quarks prefer to be bound together 
by gluon exchanges in this highly-densed thermal soup. Interestingly, further investigations 
into whether colour multiquark states in general could exist within this nuclear phase give 
positive results 

It was suggested quite a while ago in Ref. [8] that it is possible to have k < A^-baryons 
in NsusY = 4 background. In the gluon-deconfined phase, since free strings solution is 
allowed in the corresponding gravity dual theories [1], the coloured states could also exist 
in the plasma. Various possibilities of exotic multiquark states are studied in Ref. [TU]- 
|13] . Colour multiquark states in the gluon-deconfined plasma are studied in Ref. [7] where 
k > A^-baryons as well as other classes of exotic multiquark states including N + A;-baryons 
and bound state of j mesons are investigated. The phase diagram of the colour multiquarks 
nuclear phase, chiral-symmetric (xS-QGP) phase, and the vacuum nuclear phase reveals that 
colour multiquarks are thermodynamically stable in the region where the temperature is not 
too high and the density is sufficiently large (Figure 8 of Ref. [7]). 

In certain situations such as in the core of the neutron stars or other enormously densed 
astrophysical objects, exceptionally strong magnetic field is produced in addition to the high 
temperature and density. Under these fierce conditions, nuclear matters are pressed together 
so tightly that deconfinement phase transition could occur. As is shown in the phase diagram 
of Ref. [7], coloured multiquark states can exist in the intermediate range of temperature 
and sufficiently high baryon chemical potential (implying high baryon density). They are 
thermodynamically preferred over the other phases such as the vacuum and the chiral- 
symmetric deconfined phase of quark-gluon plasma (xS-QGP). It is therefore interesting to 
explore magnetic properties of the nuclear phase where coloured exotic multiquarks exist 
under these extreme situations. It is possible that certain classes of densed stars are in 
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the range of temperature and density suitable for the coloured multiquarks in the gluon- 
deconfined soup and the magnetic properties of these states thus significantly determine 
their stellar structures. 

Responses of the holographic nuclear matter to the external magnetic field have been 
intensively investigated in Ref. [ll]-[19]. It was found in Ref. [H] that the external magnetic 
field makes gluon-deconfined vacuum more stable thermodynamically than the case when 
there is no magnetic field, i.e. the transition temperature into the chiral-symmetric quark- 
gluon plasma increases with the magnetic field and saturates in the limit of an infinite field. 
Authors of Ref. [19] found a phase transition induced by the external magnetic field in the 
xS-QGP phase. This could be traced back to the nonlinearity of the DBI action used to 
describe the holographic nuclear matter. Since this transition occurs when the magnetic 
field changes from small to large strength, the Yang-Mills approximation approach [T7] is no 
longer valid and similar transition is not found without consideration of the full DBI action. 
We take the full DBI approach and investigate the magnetic responses of the multiquark 
nuclear phase with broken chiral symmetry in this article. We found that the magnetized 
multiquarks phase are always thermodynamically preferred over the magnetized vacuum. 
At a fixed density, it is also found that the multiquark states can satisfy the scale fixing 
condition up to certain critical values beyond which they would change into multiquarks 
with smaller colour charges. For higher magnetic fields, all of the multiquarks cannot satisfy 
the scale fixing condition at the same density and we would expect other phases to set in 
or the density has to be increased for the multiquark configuration to be able to satisfy the 
scale fixing condition. 

There are two multiquark configurations found below a critical field. The two configura- 
tions merge into one at the critical field and temperature for a fixed density. By comparing to 
the pure pion gradient and the xS-QGP phase, the multiquark phase is found to be preferred 
thermodynamically at large densities and moderate fields. 

In Sec. 2, the essential features of the multiquarks are reviewed. Magnetic responses 
and relevant magnetic phases of the colour multiquarks are studied in Sec. 3 using the DBI 
action. Comparison to the pure pion gradient and the xS-QGP phase is discussed in Sec. 4. 
We discuss the results and make some conclusions in Sec. 5. 

2 Exotic multiquark states in the Sakai-Sugimoto model 

In the Sakai-Sugimoto model, gluon deconfinement and the chiral-symmetry restoration are 
two distinct phase transitions. Generically they occur at different temperatures. When the 
gluons become deconfined at the deconfinement phase transition, quarks could still be bound 
together by the free gluons due to the fact that the couphng is still strong (provided that the 
density is sufficiently high) and therefore the chiral symmetry could still be broken. Due to 
the deconfinement, the bound states of multiquarks are not colour singlet in general. Certain 
properties of the coloured multiquarks are studied in Ref. [7] where it is demonstrated that 
the coloured states could exist with thermodynamical stability. When the temperature rises 
further, the bound states become less and less stable and finally completely dissolved into the 
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quark-gluon plasma. The chiral symmetry is restored and everything becomes completely 
deconfined. 

It was proposed by Witten [SOj, Gross and Ooguri [21] that a D-brane wrapping internal 
subspace of a holographic background could describe a colour-singlet bound state of quarks 
in the dual U{N) gauge theory. A wrapping D-brane sources U{1) gauge field on its world 
volume and induces an N units of f/(l) charge upon itself. This charge needs to be cancelled 
by external strings connecting to the wrapping brane. The wrapping brane with A^ strings 
attached is called a baryon vertex. 

In the gluon-deconfined phase, more strings can be attached to the baryon vertex provided 
that there are equal number of strings stretching out and go to the background horizon. This 
configuration still conserves the U{1) charge of the brane and solve the equation of motion 
of the Nambu-Goto action [H]. We can parameterize the number of radial strings stretching 
from the vertex to the horizon as kr and the number of strings connecting the vertex to 
the boundary of the background as kh. For the k > A^-baryon, kh — kr = N whilst for 
k < A^-baryon, kh + kr = N. Other classes of exotic multiquark states can be constructed 
by adding more strings in and out of the vertex. Few examples are given in Ref. [7] where 
some interesting properties are also discussed. 

There could exist an interaction among the multiquarks in the form of connecting strings 
between each vertex very similar to the string connecting two end points of quark and anti- 
quark in the holographic meson configuration. A multiquark can use one of the radial strings 
to merge with another radial string from neighbouring multiquark and form a colour binding 
interaction (while keeping kh fixed). Therefore the number of radial strings represents the 
colour charges of the multiquark. When the gluons are deconfined, the "direct" colour inter- 
action would be approximately the same as the meson and baryon binding potential of the 
Coulomb type plus some screening effect. Neglecting the direct interaction and considering 
only the DBI- induced collective behaviour of the gas of multiquarks [6l [7], an approximate 
phase diagram can be obtained showing exotic nuclear phase where multiquarks can ex- 
ist with thermodynamic stability. Schematic configurations of the three gluon-deconfined 
phases are given in Fig. [1] where the direction along the circle is the compactified coordinate 
X4 and the vertical direction is the radial coordinate u. 

3 Magnetic properties of the coloured multiquarks in 
the nuclear phase 

The setup we use is the Sakai-Sugimoto (SS) model with the source terms from the instanton 
embedded in the D8 — D8 configuration, and the radial strings similar to the configuration 
used in Ref. [7]. The instanton (the baryon vertex being pulled up all the way to the position 
of the D8-branes by the strings connecting between the vertex and the fiavour branes) is 
embedded within the D8-branes and acts as a source for the baryon density, d. The radial 
strings stretching from the instanton down to the horizon of the background act as another 
source. The number of radial strings is parameterized by = (number of radial strings) /A^. 
It also tells us how much colour charges a multiquark has. 
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(a) (b) (c) 



Figure 1: Configurations of xS-QGP (separate D8, D8)(a), vacuum (merging D8 and D8)(b) 
and exotic nuclear pliase (vertex attaclied to tlie D8-D8 witli radial strings stretcli down to 
liorizon)(c). 



The baryon chemical potential is also generated on the D8-branes by the vector part, , 
of the U{1) subgroup of the U{nf) flavour group of the D8-branes. The magnetic field is 
then turned on by another part of the U{1). We choose the direction of the magnetic field 
so that the vector potential is 



The vector part is related to the baryon chemical potential by 

/i = a^iu — oo), 

1 F)'^ 

'^'-'source /cf\ 

fisource = ■ [^) 

The contributions from the sources, fisource, are from the baryon vertex and the radial strings. 
The full expressions are given in the Appendix. The contribution from the U (1) vector gauge 
field in the D8-branes, fi, corresponds to the baryon chemical potential from the content of 
the plasma. The five- dimensional Chern-Simon term of the D8-branes generates another 
axial part of the U{1), a^, by coupling it with B and Uq . In this way, the external magnetic 
field induces the axial current Ja associated with the axial field a^. 
The background metric of the Sakai-Sugimoto model is 

ds' = (^-^y {f{u)dt' + S.jdx'dx^ + dx,') + (u'dnl + 



2'kN ^ ( u 



3/4 



i^(4) = -r^e4, , Rm = ^QsNlt, 

1/4 \riDiJ 

where f{u) = 1 — u^/u^, uj- = IGyr^i^p^T^/O. The volume of the unit four-sphere ^4 is 
denoted by V4 and the corresponding volume 4-form by £4. Is and Qs are the string length 
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scale and the string coupling. The 0:4 coordinate is compactified with radius R which is 
generically different from the curvature Rd^^ of the background. 

The DBI and the Chern-Simon actions of the D8-branes in this background can be 
computed to be 



Scs 



Af 



I I 



2 ./,, _ 



(3) 
(4) 



The normalization factor, M = NR'j~,^/{Q7r'^(27ra')^), represents the brane tension. The 
explanation of the factor 3/2 is given in Ref. [16] where it could be understood as representing 
the edge effect of the finite region with uniform magnetic field. The equations of motion 
with respect to Qq , af are 
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v/l + /(w)(af)2-(a; 
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(5) 
(6) 



The quantities d,jA are the corresponding density and current density at the boundary of 
the background {u — )• 00), they are defined to be 



j^[x, u — 00 J 



(7) 



Explicitly, they are 
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3a 



^w> + B 
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^l + f(u ){a'fy-(a 'vy + f(u)u^: 



(9) 
(10) 



For our multiquark configuration, the D8-branes starts from u = Uc and extends to m — > 00. 
At the boundary (n — )■ 00), the chiral symmetry is broken and therefore the value of af{oo) 
is taken to be a physical field, \/ip [16], describing the degree of chiral symmetry breaking. 
The total action is minimized with respect to af{oo) if the axial current Ja (also defined at 
the boundary) is zero. 

The total action does not depend on X4^{u) explicitly, therefore the constant of motion 
leads to 



{xMf 



u'[fiu)iCiu) + DjuY) - iu - jBfi + 3 Ba^, 
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where 



Wc/7K)a//K)(CK) + D{u,y) - {jA - f 5/i + 3Eao^K))2 x',{u,) 

F = ^- (12) 

v/1 + f{u,)ul x'iiu,) 

and C{u) = + B'^u'^, D{u) = d + 3Baf{u) — 3B\/(p/2. The expression of x'^^i^Uc) is given 
in the Appendix. It is determined from the force condition and the scale fixing condition 

POO 

Lo = 2 x'^{u) du = 1. (13) 



Since x'^{u) depends on both aj'(u), a^{u), we need to solve the differential equations ([5]) 
and iQ with x'^{u) substituted into the equations of motion and check whether the solutions 
satisfy the scale fixing condition Eqn. ( |T3l) . The values of the vector and axial field at 
the vertex are also chosen so that a^(Mc) = fJ'sourccttiiuc) = 0. We basically perform the 
shooting algorithm by choosing the value of /i and VV^ the expression for x'^^luc) until we 
hit a^(oo) = fi and a^(oo) = \/(p. If the resulted solution satisfies the scale fixing condition 
Lq = 1, we keep the solution. If not, we adjust the value of Uc and perform the shooting 
procedure again. The position Uc for = is given as a function of the density, the magnetic 
field, and the temperature in Fig. O 

From the solutions of the equations of motion, the relations between baryon chemical 
potential (/x) and the baryon density [d), the magnetic field [B), and the temperature (T) 
are obtained for the choice of parameters Ug = (normal baryon), = 0, as are shown in 
Figure |3l There are two types of solution corresponding to the two holographic multiquark 
configurations. One is the configuration with Uc close to Ut (configuration-A) and another 
is the configuration with a large separation between Uc and ut (configuration-B). 

The baryon chemical potential is found to be an increasing function of the density for 
most range of d for both configuration A,B. As is shown in Fig. [2l configuration-A has 
the position of vertex Uc closer to the horizon ut than configuration-B. At very small d, 
the two configurations emerge separately as two distinct configurations. Interestingly, as 
the magnetic field and temperature increase, the two configurations converge into a single 
configuration as we can see the position Uc approaches the same value at the critical field and 
temperature (see Fig. [2]). However, when the two configurations merge, the configuration no 
longer satisfies the scale fixing condition Lq = 1 and we expect it to change into other phases 
such as the chiral-symmetric quark-gluon plasma for a fixed density. It turns out that if the 
density is allowed to change, the multiquark configuration can continue to satisfy the scale 
fixing condition at higher fields provided that the density is sufficiently large. This will be 
discussed more in Section HI 

In Fig. [3l the baryon chemical potential is an increasing function of d, this is true for 
both configuration-A and B. It is roughly a linear function of the density, showing that the 
DBI-induced collective interaction between the multiquarks are negligible. As d gets larger, 
the DBI-induced effect sets in and the negative binding interaction makes /j, grows with d 
less quickly than the linear progression. Note that this DBI-induced interaction occurs even 
when the baryon is colour singlet due to the nonlinear nature of the DBI action. The origin 
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Figure 2: Position of the vertex for Ug = (normal baryon) and fixed Ja = as a function 
of (a) d with fixed B = 0.10, T = 0.10,(b) B with fixed d = 1,T = 0.10,(c) T with fixed 
B = 0.10, d = 1. The lower (blue) line is the configuration-A with u^. close to and the 
upper (red) line is the configuration-B with large separation between and ut- 
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Figure 3: The baryon chemical potential /i for = (normal baryon) and fixed ja = as a 
function of (a) d with fixed B = 0.10, T = 0.10, (b) B with fixed d = l,T = 0.10,(c) T with 
fixed B = 0.10,d = 1. The upper (blue) line is the configuration- A with Uc close to ut and 
the lower (red) line is the configuration-B with large separation between Uc and ut- 
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Figure 4: The gradient of the scalar field \jip = a^(oo) for = (normal baryon) and fixed 
jA = as a function of (a) d with fixed B = 0.10, T = 0.10,(b) B with fixed d = 1,T = 
0.10, (c) with fixed B = 0.10, T = 0.10, (d) T with fixed B = 0.10, d = 1. The lower (blue) 
line is the configuration- A with Uc close to ut and the upper (red) line is the configuration-B 
with large separation between and ut- 
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of this DBI-induced interaction is the tidal weighf of the DBI action contributed from both 
the branes' worldsheet metric and the background gauge field strength. Naturally, any form 
of energy contributes to the tidal weight even the colour singlets. 

For configuration-B, there seems to be minimal density dmin below which the shooting 
algorithm could not find other valid solutions. We are not certain what happens below these 
values. It is possible that when the field is turned on, the D8-branes acquire higher tension 
and therefore the configuration requires minimal density to pull it down in order for the 
distance between D8 and D8 to reach Lq = 1. For T = 0.10,5 = 0.10, = 0, the value of 
dmin for multiquark configuration B is approximately 0.086. 

Figure [3] shows that the chemical potential is a decreasing function with respect to the 
magnetic field. This is similar to the behaviour of baryons in chiral-symmetric quark-gluon 
plasma studied in Ref. (TB]. When the field gets stronger up to certain values, the field 
becomes too strong for the force condition to hold at the scale fixing Lq = 1. This strange 
behaviour is shown in Fig. [6] where multiquarks with smaller Ug are shown to be able to exist 
up to stronger fields. 

As is also shown in Fig. El the relationship between fi and T is as we expect, a decreasing 
function of T for fixed density d since higher temperature will melt the multiquarks away. 
For fixed d and B, the multiquark configuration satisfies the scale fixing condition up to a 
maximum temperature above which we expect it to melt into the plasma. For Ug = 0, this 
critical temperature is about 0.159 for (i = 1. 

It is interesting to note that the baryon chemical potential of = multiquarks for both 
configurations converge to the same value at critical field (~ 0.63) and temperature 0.159) 
for d = 1. This behaviour also shows up in the gradient scalar field as is shown in Fig. |H 

Fig. m shows the relations between the field VV^ the density, the magnetic field, 
the baryon chemical potential, and the temperature. The pion gradient \/ip represents the 
domain wall of the scalar field induced by the magnetic field on the nuclear vacuum |22] . 
Roughly speaking, it quantifies the degree of chiral symmetry breaking. The domain wall car- 
ries baryon charge and thus contributes to the baryon density. For multiquark configuration- 
B, it increases with B for a fixed density. From Fig. HJ the pion gradient field increases 
linearly with respect to the field for small fields. Then it starts to saturate closed to the 
critical field. This is somewhat similar to the behaviour of the pion gradient in the con- 
fined phase studied in Ref. [TB]. For configuration-B, the pion gradient field is a decreasing 
function of the density when the field is fixed. This implies that for a fixed magnetic field, 
the population of the domain wall becomes lesser as the density of the baryon (including 
multiquarks and other bound states) increases. We will see this behaviour again in Sect. 4 
when we consider the pure pion gradient phase. Finally from Fig. HI the degree of chiral 
symmetry breaking VV decreases as temperature rises for multiquark configuration-B. 

For configuration- A, the pion gradient field decreases at first for small magnetic fields, 
but turns to rise with the field around i? ^ 0.16 until it converges to configuration-B at the 
critical field. The dependence of the field \/(p on the density at a fixed B = 0.10, T = 0.10 
shows a minimum at (i ~ 0.7, corresponding to /i ~ 1.18. Then as the density grows, the 
pion gradient increases and saturates, implying limited contribution of the domain wall for 
large baryon density. The temperature dependence of the pion gradient field for multiquark 
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configuration- A shows some peculiar behaviour. First, it becomes more negative at low 
temperatures then turns to rise and converge to configuration-B at the critical temperature. 

Figure [7] shows how the pion gradient field VV^ varies with the magnetic field B for 
Us = 0, 0.10, 0.20. For the same B, multiquarks with higher Ug responds more to the magnetic 
field by inducing larger X/^p, implying higher degree of chiral symmetry breaking. The pion 
gradient field for both multiquark configuration- A, B forms a butterfiy-wing shape graph for 
each Ug ■ The edge of the wing is at the critical field where the configuration converges and 
barely satisfies the scale fixing condition. 

The magnetization of the multiquarks nuclear matter can be defined using the regulated 
free energy, J^e = ^(/^, B) + fid, in the canonical ensemble as 



dB 

where fl{fi,B) = S[aQ(u),ai{u)]{e.o.m.) — S'[magnetized vacuum]. The action with g'q 
eliminated is given by S[aQ{u),ai{u)]{e.o.m.) = Sns + Scs where 



(14) 



Sd8 = M j du C{u), 

J Ur 



fiu){C{u) + D{u)^) - {jA - iBfi + SBa^y ' 



and Scs is given in the Appendix. The grand canonical potential is regulated with respect 
to the magnetized vacuum. The action of the magnetized vacuum with non-vanishing is 



5" [magnetized vacuum] = / ^J C{u){l + f{u)u^x'l 



du, 



where 



X'4^{u)\^ac = I ■ (15) 

^nm [fiuoWiduo) 



The position uq where x'^^Iuq) = 00 of the magnetized vacuum can be solved numerically 
from Lq = 1 (with uq replacing Uc in the limit of integration). The relation between uq and 
the magnetic field is shown in Fig. |5]for T = 0.10. As the magnetic field gets stronger, the 
position of the lowest position of the D8-D8 configuration, uq, becomes larger, in order to 
satisfy the condition Lq = 1 (implying heavier branes due to magnetic field energy). At 
T = 0.10, position of uq saturates to the value of about 1.23 (The number changes with 
temperature, of course) in the limit of an infinite field. 

The magnetization of the multiquark nuclear matter is shown in Fig. [8] for = (red), 
0.10 (green), 0.20 (blue). The magnetization is positive and increases as B increases until 
the field is close to the critical value then it starts to drop. Generically, configuration-A of 
multiquarks has larger magnetization than configuration-B. For the configuration-B (A), 
multiquarks with higher (lower) rig have higher magnetizations. As the magnetic field gets 
stronger beyond the critical field for each n^, the multiquarks will undergo a transition into 
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Figure 5: Relation between uq and external magnetic field B of the vacuum for the temper- 
ature T = 0.10, Uq saturates to the approximate value of 1.23 at large field. 
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Figure 6: Comparison between the baryon chemical potential as a function of B at fixed 
Ja = 0,d = 1,T = 0.10 and (a) = (normal baryon), the bottom graph, (b) = 0.10, the 
middle graph, (c) rig = 0.20, the top graph. The upper (blue) lines are the configuration-A 
with Uc close to ut and the lower (red) lines are the configuration-B with large separation 
between Uc and ut- 
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Figure 7: Comparison between the gradient of the scalar field VV^ a function of B at fixed 
= 0,d = 1,T = 0.10 and (a) = (normal baryon), the bottom graph, (b) Ug = 0.10, the 
middle graph, (c) Ug = 0.20, the top graph. The lower (blue) lines are the configuration- A 
with Uc close to ut and the upper (red) lines are the configuration-B with large separation 
between Uc and ut- 




Figure 8: The magnetization of the multiquarks nuclear matter at fixed = 0, c? = 1, and 
T = 0.10 for Ug = (red), 0.10 (green), 0.20 (blue). The upper lines are the configuration-A 
with Uc close to ut and the lower lines are the configuration-B with large separation between 
Uc and Ut- 
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Figure 9: The free energy of the multiquarks nuclear matter at fixed Ja = 0,d = 1, and 
T = 0.10 for rig = (red), 0.10 (green), 0.20 (blue). The upper lines are the configuration-A 
with Uc close to ut and the lower lines are the configuration-B with large separation between 
Uc and Ut- 

the ones with smaller Ug. For even larger fields, even the Ug = multiquarks cannot satisfy 
the scale fixing condition if the density is not allowed to change. 

Interestingly, numerical studies reveal that the grand canonical potential of the multi- 
quark phase is always lower than the grand canonical potential of the magnetized vacuum, 
i.e. S[ao{u), ai{u)]{e.o.m.) — ^[magnetized vacuum] < 0, for the entire range of B. This sug- 
gests that once /i > fionset, the magnetized multiquark phase is always thermodynamically 
preferred over the magnetized vacuum, the situation similar to the case when there is no 
magnetic field investigated in Ref. [7]. Among the two configurations, we found from Fig. M 
that the free energy of configuration-B is always lower than configuration-A and thus more 
stable thermodynamically. These two multiquark configurations- A, B are the long and short 
cusp configurations discussed in Ref. [6], being extended to the general case with nonzero 
magnetic fields. It is found here that for a fixed density, strong field and/or high tempera- 
ture (see Fig. [3] and H]) converge the two into a single configuration right before dissociating 
them altogether. 

Figure [9] shows how the free energy changes with the magnetic field for tt.^ = (red), 
0.10 (green), 0.20 (blue) at the temperature T = 0.10 and the density d = 1. For each 
Ug, both configurations converge to the same configuration (with the same baryon chemical 
potential, degree of chiral symmetry breaking and free energy) at the critical fields. The 
critical fields for Ug = 0, 0.1, 0.2 are roughly 0.63, 0.48, 0.38 respectively. 

4 Comparison to other phases 

In this section we compare the baryon chemical potential and free energy of the magne- 
tized multiquarks to the pure pion gradient phase and the chiral symmetric quark-gluon 
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plasma (x5-QGP) phase, both under the external magnetic field with gluons deconfined. 
The pure pion gradient phase is defined to be the phase with figource = (sourceless case) 
and the baryon chemical potential comes purely from the induced gradient field, \/(p, in 
response to the external field. The baryon density also comes purely from the pion gradient 
field {d = 3B'\/(p/2). A similar situation in the confined phase of the antipodal SS model has 
been studied in Ref. [IS]. The Xs-QGP under the presence of the external magnetic field has 
been explored in Ref. [121 [JH] but again only limited to the antipodal case of the SS model. 
In this section we explore some of their magnetic properties in more general case where x'^^lu) 
is not zero and the scale is fixed to Lq = 1. Even though the extra constraints are irrelevant 
to the xs-QGP (since x'^ = for this configuration), it makes crucial difference in the case 
of pure pion gradient phase. The scale fixing condition is found to be very difficult for the 
pure pion gradient configuration to satisfy for most of the density as we will discuss below. 

All three phases under consideration obey the same set of equations of motion, Eqn. ([5]), ([6]) 
with each specific set of the boundary conditions and parameters as the following, 

multiquark phase : Ja = 0,Hsource = a^{uc),\7^p = af{oo),af{uc) = 0, 

pure pion gradient phase: all the same with the multiquark phase with the following ex- 
ceptions, flsource = 0, a^(uc) ^0,d= ^Bx/if, 

Xs-QGP : x'^{u) = and \/(p = a^{oo) = O,fisource = do {u^ = Ut) = OJa = \Bii (since 
the configuration extends to ut and /{ut) = so that Eqn. ([5]) is zero). 

First, we will explore certain properties of the pure pion gradient phase and show that 
it does not exist in the range of parameters (ci > 1, -B < 1 — 2) under consideration. Then 
comparison between the multiquark and the xs-QGP phases will be discussed. 

4.1 Pure pion gradient phase 

For pure pion gradient configuration, the contribution of the sources, the vertex and strings, 
is set to zero.. Effectively, we set fisource = 0,d = 3B\/(p/2. This is because when VV^ is 
zero, the density d should represent the density of the sources, i.e. the pure multiquark or 
pure baryon configuration, therefore the source density should be given by ci — ^Bs/ip on the 
right-hand side of Eqn. When we fix the value of the density at a fixed magnetic field, 
\/ip is also fixed. For example, when d = 1, B = 0.1, s/ip ~ 6.667, a relatively large value. 
This large value of VV^ = a^(oo) leads to a generically large value of af{u) for the most 
range of u. From Eqn. ( ITT]) and (fT2]) . we see that for the pure pion phase, D{u) = ^Ba^iu) 
and thus it must be large for the most range of u as well. In the multiquark configuration, 
the d dependence of D{uc) in the expression of F, Eqn. (fT2]) . will compensate the largeness of 
D{u) and X4 can be made sufficiently large so that Lq = 1 could still be satisfied. However, 
in pure pion phase, D{uc) is simply zero. This makes X4 getting smaller as the density gets 
larger and the scale fixing condition Lq = 1 would not be satisfied above certain value of the 
density for a fixed B. 
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Figure 10: The position Uc, the pion gradient field, and the baryon chemical potential of the 
pure pion gradient phase aX d = 0.10, T = 0.10 as a function of the magnetic field. 



As a result, we wish to keep \/ip sufficiently small in order to satisfy the scale fixing 
condition. This implies that higher densities require larger magnetic fields. To demonstrate 
this, we fix baryon density to d = 0.1 and plot the position Uc of the vertex and the baryon 
chemical potential as a function of the magnetic field in Fig. [TOl The graph of Uc shows a 
minimal field at about B ^ 0.229 below which Lq < 1 for all solutions. For a larger density 
d > 1, the required field strengths are i? >> 1 in order for the scale fixing condition to be 
satisfied. For the range of parameters d = 1.0, B < 1.0, we therefore need to consider only 
the two phases of the multiquark and the X5-QGP. The same situation occurs for the range 
of parameters d = 10, B < 1 — 2 where the pure pion gradient phase does NOT satisfy the 
scale fixing condition and therefore does not exist as well. 

4.2 Multiquark-domain wall versus xs-QGP phase 

The baryon chemical potential fi is to be found by shooting algorithm for a fixed d, B, T, for 
each phase. For d = 1, B = 0.10, T = 0.10, they are shown in Fig. [TTl Observe that there 
are two possible solutions for the xs'-QGP phase. As the magnetic field increases beyond 
a certain value (in this case around B 0.25), there will be phase transition to another 
solution within this phase. This behaviour is explored in details in Ref. [12]. When the 
density is raised to d = 10, the transition occurs at higher field around B ^ 0.86 (Fig. [T^ . 
The transitions can also be seen in the plots of the free energy. Fig. [T2][T5l where the slopes 
of the graphs change abruptly around the critical fields. For d = 1, this is quite small and 
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Figure 11: Comparison between the baryon chemical potential for T = 0.10 at a fixed density 
d = 1 of (a) Ug = (normal baryon) multiquark configuration- A, the top (blue) graph, (b) 
Xs-QGP, the middle (black) graph, (c) Ug = (normal baryon) multiquark configuration-B, 
the bottom (red) graph. 




Figure 12: Comparison between the free energy for T = 0.10 at a fixed density d = 1 
of (a) Ug = (normal baryon) multiquark configuration- A, the top (blue) graph, (b) Xs- 
QGP, the middle (black) graph, (c) = (normal baryon) multiquark configuration-B, the 
bottom (red) graph . 
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Figure 13: The position of the vertex function of -B at a fixed density d = 1 (lower) 

and d = 10 (upper) for T = 0.10 of = (normal baryon) multiquark configuration-B 
phase. 
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Figure 14: Comparison between the baryon chemical potential for T = 0.10 at a fixed density 
c? = 10 of (a) Us = (normal baryon) multiquark configuration- A, the top (blue) curve, (b) 
Xs-QGP, the black curve, (c) Ug = 0.2 multiquark configuration-B, the dashed red curve, (d) 
Us = (normal baryon) multiquark configuration-B, the red curve. 
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Figure 15: Comparison between the free energy for T = 0.10 at a fixed density d = 10 
of (a) Us = (normal baryon) multiquark configuration- A, the top (blue) curve, (b) xs- 
QGP, the black curve, (c) Us = 0.2 multiquark configuration-B, the dashed red curve, (d) 
Us = (normal baryon) multiquark configuration-B, the red curve. 




Figure 16: Plots between the pion gradient field of the multiquark phase and the magnetic 
field for T = 0.10, ra^ = at d = 1 (shorter) and c? = 10 (longer). 
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somewhat hard to see but it becomes apparent for d = 10. 

From the plots of the free energy, Fig. [T21 the multiquark configuration- A is the least 
preferred phase when the density is small {d = 1). Its free energy is larger than the x^-QCP 
phase for all fields. For B < 0.196, the most preferred phase is the multiquark configuration- 
B phase with the lower free energy. The xs-QGP phase is more stable for d = 1, B > 0.196. 
Nevertheless, the multiquark configurations can exist up to only about the critical fields 
beyond which they cannot satisfy the scale fixing condition at that particular density. 

However, this does not mean that the multiquarks phase cannot exist in the range of 
field larger than the critical value. Stronger field gives the D8-branes larger tension and thus 
it requires sufficiently heavier vertex and strings to pull it down in order for the distance 
between D8 and D8 to reach Lq = 1. This implies that we need larger d in order to make 
the configuration satisfy the scale fixing condition at stronger fields. Fig. [T3|14II15I confirm 
this insight. They show the plots of the multiquarks configurations when the density is 
large {d = 10). Multiquark configurations can exist far beyond the critical field B ^ 0.63 of 
the small d case {d = 1). In particular. Fig. [15] demonstrates that at ci = 10, the multiquark 
configurations (n^ = 0,0.2), with lower free energies, are thermodynamically preferred over 
the X5-QGP for B < 0.61 and B < 0.348 respectively. 

It is thus reasonable to conclude that for larger densities, the multiquarks phase will be 
more and more preferred over the x^-QGP phase, in a larger and larger range of the field. 
Magnetized multiquarks and the induced pion gradient field are thus stable thermodynami- 
cally and they will mix together in the magnetized nuclear (multiquark-domain wall) phase 
provided that the density is sufficiently large and the temperature is not too high. 

Finally for completeness, we present the plots of the pion gradient field of the multi- 
quark phase (Fig. [T6|l . the pion gradient field becomes smaller for a given B as the density 
increases. However, it extends to larger range of fields for larger density. We can therefore 
conclude that at the large densities (and baryon chemical potential), contribution of the pion 
gradient becomes lesser and the multiquarks contribute dominantly to the baryon density 
and chemical potential. This is also shown in Fig. m 

5 Discussions and Conclusion 

In Sakai-Sugimoto model, chiral symmetry restoration and gluon deconfinement are two 
distinct phase transitions. Generically, with an exception of the antipodal case with X4 = 0, 
gluon deconfinement occurs at lower temperature than the chiral symmetry restoration. For 
the region of the phase diagram between the two transitions, coloured multiquarks can exist 
with thermodynamical stability (the phase diagram is shown in Figure 8 of Ref. [7]). 

Magnetic responses of the nuclear phase with colour multiquarks are studied here by 
using one component of the U{1) subgroup of U{Nf) as the vector potential of the external 
magnetic field. The Chern-Simon action of the D8-branes couples the magnetic field to 
an axial vector component, a^, of the U{1), inducing axial current Ja- When the chiral 
symmetry is broken, we effectively set Ja to zero. The value of af(oo) then describes the 
degree of chiral symmetry breaking of the phase. 
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There are two possible multiquark configurations A and B. Configuration- A is the config- 
uration where the baryon vertex is close to the horizon. Configuration-B, on the other hand, 
has the baryon vertex more separated from the horizon. By comparing the free energy of the 
two configurations in Fig. |9l we found that configuration-B is more stable themodynamically. 
We establish relations between the baryon chemical potential and the baryon density, the 
external magnetic field, and the temperature for both configurations as are shown in Fig. Ei 
Baryon chemical potential is an increasing function of the density when the field is turned 
on. This is the same behaviour to the case when there is no field. 

On the other hand, the relation between chemical potential and the magnetic field is 
rather interesting. The baryon chemical potential is a decreasing function of the field. For 
multiquarks with high value of (number of radial strings), the configuration finds it more 
difficult to satisfy the scale fixing condition at large fields. There is a maximum field strength 
for each Us above which the multiquark configuration cannot exist (Fig. [6]). This is in contrast 
to the behaviour of the chiral-symmetric quark-gluon plasma (in the antipodal case of the 
Sakai-Sugimoto model with no instantons, i.e. x'^{u) = case) studied in Ref. [16] where 
chemical potential is always a decreasing function with respect to B and the configuration 
continues to exist at arbitrarily large fields. This is due to the fixation of the density. Stronger 
field gives the flavour branes more tension and when the field is too strong, a fixed density 
source would not be sufficiently heavy to pull the branes down for the distance between D8 
and D8 to reach Lq = 1. Temperature also has effect on the multiquarks, sufficiently high 
temperature will melt away the multiquarks even in the presence of an external field. 

The gradient of the scalar field representing the chiral symmetry breaking, \/(p = a^(oo), 
is found to roughly increase in magnitude with the field. For the same field strength and 
fixed density, multiquarks with higher Ug (i.e. larger colour charges) show higher degree 
of chiral symmetry breaking (larger magnitude of a^(oo)), but can only sustain the force 
condition up to smaller fields as is shown in Fig. [71 

The mixing of pion gradient with the miltiquark in the multiquark phase decreases as 
the density increases (Fig. Hj). It is found that the pure pion gradient phase (no multiquark 
contribution) does not satisfy the scale fixing condition for large densities and moderate 
fields. 

What would happen if the magnetic field increases beyond the point where the multi- 
quarks can satisfy the scale fixing condition = 1? We would expect the multiquarks to 
change into the multiquarks with lower as is shown in Fig. O for a fixed d and T since 
they can still satisfy the scale fixing condition. This induces a sudden drop in the baryon 
chemical potential. Also in the situation where /i is kept fixed instead of d, the multiquarks 
are forced to jump to the larger d in order to change into the multiquarks with lower Ug as 
the field increases beyond the critical point. For even larger fields, all of the multiquarks 
cannot satisfy the scale fixing condition for a fixed density. There would be phase transition 
to other phase. For a fixed density, the phase will change into the X5-QGP. However, if 
we allow the density to change (in a more realistic situation), the system could change into 
the multiquark (with pion gradient mixing) phase for a sufficiently large density. The phase 
of multiquark with pion gradient mixing is found to be more preferred than the Xs-QGP 
at large densities (implying large baryon chemical potentials) and moderate fields. This is 
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shown in Fig. [13 

For configuration-B multiquarks, The magnetization of the multiquark nuclear matter 
is found to be an increasing function of B for Ug = 0,0.10,0.20 except when the fields 
get close to the critical points. Close to the critical fields, the magnetizations saturate and 
even start to decrease. The magnetized multiquarks phases are thermodynamically preferred 
over the magnetized vacuum once the baryon chemical potential is higher than the onset 
value [n > fionset)- This is similar to the case when there is no magnetic field investigated 
in Ref. [7]. 
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A Force condition of the multiquark configuration 

The forces on the D4-brane in the fiavour D8-branes are balanced among three forces from 
the tidal weight of the D4-brane, the force from the strings attached to the D4, and the 
force from the D8-branes. Varying the total action with respect to Uc gives the surface term. 
Together with the scale-fixing condition 2 dux'^{u) = Lq = 1, we obtain [B] 



/ dS 



duc J / dx'^ 



(16) 



as the condition on Uc- 

The Legendre transformed action is given by 



S 




u^x'l 



X ^^\u){C{u) + D{uf)-{3A-\B^^ + Wal)\ (17) 

where C{u) = + B'^v?,D{u) = d + 3Baf{u) — 3B\/(p/2. It is calculated by performing 
Legendre transformation with respect to a^' and af' respectively. Note that the Chern-Simon 
action is also included in the total action during the transformations. 
The Chern-Simon term with the derivatives a^', a^' eliminated is 

3 .oo (a^ijA - iBfi + 3Ba^) - f{u)D{u)af) Jjh+u^xf 
Scs = ~M-B I du ^ " 



f{u){C{u) + D{uy) - (^jA - iBfi + 3Bal 
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Lastly, in order to compute x'^{uc) we consider the source term [7| 



+ ns{uc - ut) 



where Ug = k^/N is the number of radial strings in the unit of 

From Eqn. ([16]) , ([IT]) , (HI]) , ([20]) , and setting al{u,) = source. a^{u,) 
obtain 



(19) 
(20) 

we can solve to 



(4K))' 



fcU. 



9 {aC, + Dl)~{jA-\B^i 



rf2 



- 1 



where fc = f{uc),Cc = C{uc),Dc = D{uc). 

When we fix the parameter n^, the temperature T, the baryon density d, the axial current 
jA = (by minimizing the action with respect to af(oo)), and setting af(Mc) = 0,a^(Mc) = 
fJ'source, then the position Uc of the D4-brane is completely determined as a function of the 
magnetic field B. Once the equations of motion are solved, the value of /i = aj'(oo) and 
a{^(oo) are determined. 
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